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Abstract 

In this paper, we introduce a multivariate extension of the classical univariate Value-at-Risk 
(VaR). This extension may be useful to understand how solvency capital requirement is affected 
by the presence of risks that cannot be diversify away. This is typically the case for a network 
of highly interconnected financial institutions in a macro-prudential regulatory system. We also 
generalize the bivariate Conditional- Tail-Expectation (CTE), previously introduced by Di Bern- 
ardino et al. (2011), in a multivariate setting and we study its behavior. Several properties 
have been derived. In particular, we show that these two risk measures both satisfy the positive 
homogeneity and the translation invariance property. Comparison between univariate risk meas- 
ures and components of multivariate VaR and CTE are provided. We also analyze how they are 
impacted by a change in marginal distributions, by a change in dependence structure and by a 
change in risk level. Interestingly, these results turn to be consistent with existing properties on 
univariate risk measures. Illustrations are given in the class of Archimedean copulas. 

Keywords: Multivariate risk measures, Level sets of distribution functions, Multivariate 
probability integral transformation, Stochastic orders, Copulas and dependence. 



Introduction 

During the last decades, researchers joined efforts to properly compare, quantify and manage 
risk. Regulators edict rules for bankers and insurers to improve their risk management and to 
avoid crises, not always successfully as illustrated by recent events. 

Traditionally, risk measures are thought of as mappings from a set of real-valued random vari- 
ables to the real numbers. However, it is often insufficient to consider a single real measure to 
quantify risks created by business activities, especially if the latter are affected by other external 
risk factors. Let us consider for instance the problem of solvency capital allocation for financial 
institutions with multi-branch businesses confronted to risks with specific characteristics. Under 
Basel II and Solvency II, a bottom-up approach is used to estimate a "top-level" solvency capital. 
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This is done by using risk aggregation techniques who may capture risk mitigation or risk diver- 
sification effects. Then this globai capital amount is re-allocated to each subsidiaries or activities 
for internal risk management purpose ("top-down approach"). Note that the solvability of each 
individual branch may strongly be affected by the degree of dependence amongst all branches. 
As a result, the capital allocated to each branch has to be computed in a multivariate setting 
where both marginal effects and dependence between risks should be captured. In this respect, 
the "Euler approach" (e.g., see Tasche, 2008) involving vector- valued risk measures has already 
been tested by risk management teams of some financial institutions. 

Whereas the previous risk allocation problem only involves internal risks associated with busi- 
nesses in different subsidiaries, the solvability of financial institutions could also be affected by 
external risks whose sources cannot be controlled. These risks may also be strongly heterogeneous 
in nature and difficult to diversify away. One can think for instance of systemic risk or contagion 
effects in a strongly interconnected system of financial companies. As we experienced during the 
2007-2009 crisis, the risks undertaken by some particular institutions may have significant impact 
on the solvability of the others. In this regard, micro-prudential regulation has been criticized 
because of its failure to limit the systemic risk within the system. This question has been dealt 
with recently by among others, Gauthier et al. (2010) and Zhou (2010) who highlights the benefit 
of a "macro-prudential" approach as an alternative solution to the existing "micro-prudential" 
one (Basel II) which does not take into account interactions between financial institutions. 

In this paper, we introduce a multivariate extension of the classical univariate Value-at-Risk 
(VaR). We also study the multivariate extension of the Conditional- Tail-Expectation (CTE), 
previously introduced by Di Bernardino et al. (2011) in a bivariate setting. These extensions 
may be useful to understand how solvency capital requirement (SCR) could be computed in a 
macro-prudential regulatory framework in which institutions can also be affected by risks under- 
taken by its competitors. Several properties have been derived. In particular, we show that these 
two risk measures both satisfy the positive homogeneity and the translation invariance property. 
Comparison results between univariate risk measures and components of multivariate risk meas- 
ures are provided. We also analyze how they are impacted by a change in marginal distributions, 
by a change in dependence structure and by a change in risk level. Interestingly, these results turn 
to be consistent with existing properties on univariate risk measures. In particular, we prove that, 
for different Archimedean families of copulas, an increase in dependence amongst risks tends to 
lower the amount of solvency capital required for these risks. At the extreme case where risks are 
perfectly dependent (i.e., comonotonic) , each component of the multivariate VaR (resp. CTE) 
are equal to the corresponding univariate VaR (resp. CTE), which in turn can be viewed as 
a lower bond. This feature may seem counterintuitive at a first sight. However, in a financial 
network where each institution bears risks that cannot be diversify away by the others, the higher 
the degree of dependence is, the smaller the number of independent (undiversified) risks is, and 
the lower the solvency capital shall be required for the global system. At the extreme case where 
risks in the financial network are perfectly comonotonic, only one single factor governs risks in 
the global system and a micro-prudential rule may be suitable for computing solvency capital. 
This is in line with the observation made by Zhou (2010): "When regulating a system consist- 
ing of similar institutions, or in other words, the system is highly interconnected, considering a 
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micro-prudential regulation can be sufficient for reducing the overall systemic risk." (Zhou, 2010). 

Furthermore, our multivariate risk measures can be considered as a "fair" allocation of solvency 
capital with respect to individual risk-taking behavior in an economy with d interconnected fin- 
ancial institutions. We prove that capital required for an institution is affected both by its own 
marginal risk and by the degree to which its business is connected to the activity of the other 
institutions. However, the solvency capital required for one particular institution does not depend 
on the marginal risks bearing by the others. This can be considered as an invariance property 
with respect to a change in external risk's marginal distributions, as far as the dependence struc- 
ture is being fixed. 

From the years 2000 onward, much research has been devoted to risk measures and many exten- 
sion to multidimensional settings have been suggested (see, e.g., Jouini et al, 2004 and Bentahar, 
2006, Henry et al, 2012). Unsurprisingly, the main difficulty regarding multivariate generaliza- 
tions of quantile-based risk measures (as the VaR and the CTE) is the fact that vector preorders 
are, in general, partial preorders. Then, what can be considered in a context of multidimensional 
portfolios as the analogous of a "worst case" scenario and a related "tail distribution"? This is 
the first question we shall address by suggesting a suitable definition of quantiles for multi-risk 
portfolios. 

In the last decade, several attempts have been made to provide a multidimensional generaliz- 
ation of the univariate quantile function. For example, Masse and Theodorescu (1994) defined 
multivariate quantiles as half-planes and Koltchinskii (1997) provided a general treatment of 
multivariate quantiles as inversions of mappings. Another approach is to use geometric quantiles 
(see, for example, Chaouch et al, 2009). Along with the geometric quantile, the notion of depth 
function has been developed in recent years to characterize the quantile of multidimensional dis- 
tribution functions (for further details see, for instance, Chauvigny et al, 2011). We refer to 
Serfling (2002) for a large review on multivariate quantiles. When it turns to generalize the 
Value-at-Risk measure, Embrechts and Puccetti (2006), Nappo and Spizzichino (2009), Prekopa 
(2012) use the notion of quantile curve (formally introduced in Section]!]). Contrarily to the latter 
approach, the multivariate Value-at-Risk proposed in this paper quantifies multivariate risks in 
a more parsimonious and synthetic way. This feature can be relevant from an operational point 
of view. 

In the literature, several generalizations of the classical univariate CTE have been proposed, 
mainly using as conditioning events the total risk or some univariate extreme risk in the portfolio. 
This kind of measures are suitable to model problems of capital allocation in a portfolio of 
dependent risks. More precisely, let X = {X\,...,Xj) be a risk vector, where, for any i = 
l,...,d, the component JQ denotes risk (usually claim or loss) associated with subportfolio i. 
Then, S = X\ + ■ ■ ■ + corresponds to the total risk of this portfolio, = min{Ai, Xj} 
and Xf<£\ = max{A"i, X^} are the extreme risks. In capital allocation problems, we are not 
only interested in the "stand-alone" risk measures CTE a (Aj) = E[Aj|Aj > Qxi( a )]> where 
Qxi(ot) = inf{x G R+ : Fx^x) > a} is the univariate quantile function of Xi at risk level a, but 
also in 
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CTE s ™(X i )=E[X i \S>Q s (a)}, (1) 

CTEr {X t ) =E[Xi | X {1) > Q Xm (a)], (2) 

CTE- ax (X 4 ) = E[Xi | X {n) > Q X(d) (a)], (3) 

for i = l,...,d. The interested reader is referred to Cai and Li (2005) for further details. For 
explicit formulas of CTE^ um (Aj) in the case of Fairlie-Gumbel-Morgenstern family of copulas, see 
Barges et al. (2009). Landsman and Valdez obtain an explicit formula for CTE™ m (JQ) in the 
case of elliptic distribution functions (see Landsman and Valdez, 2003); Cai and Li in the case 
of phase-type distributions (see Cai and Li, 2005). Furthermore, we recall that CTE™ m (Aj) is a 
key tool to calculate the amount of solvency capital in the "Euler approach" (see, e.g., Tasche, 
2008). 

The multivariate version of Conditional-Tail-Expectation, studied in this paper, is essentially 
based on a 11 distributional approach". It is constructed as the conditional expectation of a 
multivariate random vector given that the latter is located in a particular set corresponding to the 
a-upper level set of the associated multivariate distribution function (in a bivariate setting see Di 
Bernardino et al, 2011). More precisely we consider a non-negative multivariate random vector 
X = (Xi, . . . ,Xd) and the associated a-upper level set, i.e., L(a) = {x G : Fx(x) > a}, for 
a € (0, 1), where Fx. is the multivariate distribution function of X. As the total information 
of the vector X is completely described by its multivariate distribution function, the set L(a) 
captures the information coming both from the marginal distributions and from the multivariate 
dependence structure. Then contrarily to existing generalizations of the univariate CTE presen- 
ted above, our CTE Q (X) does not use an arbitrary real- valued aggregate transformation (sum, 
min, max,...). Indeed, using an aggregate procedure between the risks can be inappropriate 
to measure risks with heterogeneous characteristics especially in an external risks problem. 
Moreover, as opposed to our proposal of multivariate risk measures, the multivariate risk 
measures defined in Q-Q do not satisfy the invariance property with respect to a change in 
external risk's marginal distributions. This means that for a given risk component i in the 
portfolio, changes in marginal distributions of the other risks lead to changes in risk measures 
Q-© for name i, even if the dependence structure does not change. 

The paper is organized as follows. In Section [TJ we introduce some notations, tools and technical 
assumptions. In Section [2j we propose a new multivariate generalization of Value- at- Risk. In 
Section [3j we generalize in a multivariate setting the Conditional- Tail- Expectation, previously 
introduced by Di Bernardino et al. (2011) in dimension two. We study the properties of our 
multivariate VaR and CTE in terms of Artzner et al. (1999) 's invariance properties of risk 



measures (see Sections 2.2 and 3.1). We also compare the components of these multivariate risk 



measure with the associated univariate risk measures (see Sections 2.3 and 3.2). The behavior 



of VaR a (X) (resp. CTEq,(X)) with respect to a change in marginal distributions, a change in 



dependence structure and a change in risk level a is discussed in Sections 2.4 2.6 (resp. Sections 



3.3 3.5). Further illustrations in some Archimedean copula cases, are presented both for the 
multivariate Value-at-Risk and for the Conditional-Tail-Expectation. In the conclusion, we discuss 
open problems and possible directions for future work. 
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1. Basic notions and preliminaries 

In this section, we first introduce some notation and tools which will be used later on. 
Stochastic orders 

From now on, let Qx(ct) be the univariate quantile function of a risk X at level a € (0,1). 
More precisely, given an univariate continuous and strictly monotonic loss distribution function 
Fx, Qx(a) = F^ 1 (a), \/ a G (0, 1). We recall here the definition and some properties of useful 
univariate and multivariate stochastic orders. 

Definition 1.1 (Stochastic dominance order) Let X and Y be two random variables. Then 
X is said to be smaller than Y in stochastic dominance, denoted as X ^ s t Y , if the inequality 
Qx{ot) < Qy(a) is satisfied for all a € (0, 1). 

Definition 1.2 (Stop- loss order) Let X and Y be two random variables. Then X is said 
to be smaller than Y in the stop-loss order, denoted as X Y, if for all t € R, 

E[(X - t)+] < E[(y - *)+], with x + := max{x, 0}. 

Definition 1.3 (Increasing convex order) Let X and Y be two random variables. Then X 
is said to be smaller than Y in the increasing convex order, denoted as X <i cx Y , if E[/(X)] < 
E[/(Y")], for all non- decreasing convex function f such that the expectations exist. 

The stop- loss order and the increasing convex order are equivalent (see Theorem 1.5.7 in Midler 
and Stoyan, 2001). Note that stochastic dominance order implies stop-loss order. For more 
details about stop-loss order we refer the interested reader to Miiller (1997). Moreover, a sufficient 
condition for the stop-loss order is the dangerousness order relation as stated in the following 
lemma. 

Lemma 1.1 (Ohlin, 1969) Let X and Y be random variables with finite means such that 
K[X] < E[y], and there exists some real number c such that Fx(x) < Fy(x), for all x < c and 
Fx(x) > Fy(x), for all x > c. Then X precedes Y in dangerousness order, written X ■<£> Y, 
and this implies the stop-loss order X < s i Y . 

According to Biihlmann's terminology, when X <dY the distribution function Fy is said to be 
more dangerous than Fx ■ This terminology is essentially related to the variability of the random 
variables X and Y (see Section 3.4.2.2 in Denuit et al., 2005). For further details, the reader is 
referred to Biihlmann et al. (1977). Finally, we introduce the definition of supermodular function 
and supermodular order for multivariate random vectors. 

Definition 1.4 (Supermodular function) A function f : M rf — > R is said to be supermodular 
if for any x,y£M> d it satisfies 

f(x) + f(y)<f(xAy) + f(xVy), 
where the operators A and V denote coordinatewise minimum and maximum respectively. 
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Definition 1.5 (Supermodular order) Let X and Y be two d— dimensional random vectors 
such that E[/(.X)] < E[/( Y)], for all supermodular functions f : R" — >• R, provided the expectation 
exist. Then X is said to be smaller than Y with respect to the supermodular order (denoted by 

X Y). 

This will be a key tool to analyze the impact of dependence on our multivariate risk measures. 
Kendall distribution function 

Let X = (X±, . . . ,Xd) be a d— dimensional random vector, d > 2. As we will see later on, our 
study of multivariate risk measures strongly relies on the key concept of Kendall distribution func- 
tion (or multivariate probability integral transformation), that is, the distribution function of the 
random variable F(X), where F is the multivariate distribution of random vector X. From now 
on, the Kendall distribution will be denoted by K, so that K(a) = P[i ? (X) < a], for a G [0, 1]. 
We also denote by K(a) the survival distribution function of F(X), i.e., K(a) = P[F(X) > a]. 
For more details on the multivariate probability integral transformation, the interested reader is 
referred to Caperaa et al., (1997), Genest and Rivest (2001), Nelsen et al. (2003), Genest and 
Boies (2003), Genest et al. (2006) and Belzunce et al. (2007). 

In contrast to the univariate case, it is not generally true that the distribution function K of F(X.) 
is uniform on [0, 1], even when F is continuous. Note also that it is not possible to characterize 
the joint distribution F or reconstruct it from the knowledge of K alone, since the latter does not 
contain any information about the marginal distributions F Xl , • • ■ , Fx d (see Genest and Rivest, 
2001). Indeed, as a consequence of Sklar's Theorem, the Kendall distribution only depends on 
the dependence structure or the copula function C associated with X (see Sklar, 1959). Thus, we 
also have K(a) = P[C(U) < a] where U = (U u . . . , U d ) and C7i = F Xl (Xi), ...,U d = F Xd {X d ). 

Furthermore: 

• For a d— dimensional random vector X = (Xi, . . . , X d ) with copula C, the Kendall dis- 
tribution function K(a) is linked to the Kendall's tau correlation coefficient via: tc = 

2 2rf C i U i 1 ' f° r ^ — ^ ( see Section 5 in Genest and Rivest, 2001). 

• The Kendall distribution can be obtain explicitly in the case of multivariate Archimedean 
copulas with generatoi]^] <p, i.e., C(u\,...,Ud) = (f)" 1 (4>(u\) + • • • + 4>(ud)) for all 
(ui, . . . , Ud) 6 [0, l] d . Table [l] provides the expression of Kendall distributions associ- 
ated with Archimedean, independent and comonotonic d— dimensional random vectors (see 
Barbe et al, 1996). Note that the Kendall distribution is uniform for comonotonic random 
vectors. 

For further details the interested reader is referred to Section 2 in Barbe et al. (1996) and 
Section 5 in Genest and Rivest (2001). For instance, in the bivariate case, the Kendall 



3 Note that (f> generates a d— dimensional Archimedean copula if and only if its inverse 1 is a d— monotone on 
[0,oo) (see Theorem 2.2 in McNeil and Neslehova, 2009). 
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Copula 


Kendall distribution K(a) 


Archimedean case 

Independent case 
Comonotonic case 


«+«Eti( ln T )l ) 



Table 1: Kendall distribution in some classical d— dimensional dependence structure. 

distribution function is equal to a — , a S (0, 1), for Archimedean copulas with differ- 
entiable generator (p. It is equal to a (1 — ln(a)) , a E (0, 1) for the bivariate independence 
copula and to 1 for the counter-monotonic bivariate copula. 

• It holds that a < K{a) < 1, for all a £ (0, 1), i.e., the graph of the Kendall distribution 
function is above the first diagonal (see Section 5 in Genest and Rivest, 2001). This is 
equivalent to state that, for any random vector U with copula function C and uniform 
marginals, C(U) ^ s t C C (U°) where U c = (J7£, . . . , US) is a comonotonic random vector 
with copula function C c and uniform marginals. 

This last property suggests that when the level of dependence between X±, . . . , Xj increases, 
the Kendall distribution also increases in some sense. The following result, using definitions of 
stochastic orders described above, investigates rigorously this intuition. 

Proposition 1.1 Let U = (Ui, . . . , U^) (resp. U* = (U*, . . . ,U$)) be a random vector with 
copula C (resp. C* ) and uniform marginals. 

If U^ sm ir, then C{U)< sl C*{lT). 

Proof: Trivially, U < sm IT =>■ C(u) < C*(u), for all u G [0, l] d (see Section 6.3.3 in Denuit 
et ai, 2005). Let / : [0,1] — > M. be a non-decreasing and convex function. It holds that 
/(C(u)) < /(C*(u)), for all u G [0, l] d , and E[/(C(U))] < E[/(C*(U))]. Remark that since C* 
is non-decreasing and supermodular and / is non-decreasing and convex then / o C* is a non- 
decreasing and supermodular function (see Theorem 3.9.3 in Muller and Stoyan, 2001). Then, 
by assumptions, E[/(C(U))] < E[/(C*(U))] < E[/(C*(U*))]. This implies C(U) C*{U*). 
Hence the result. □ 

From Proposition |1.1| we remark that U ^ sm U* implies an ordering relation between cor- 
responding Kendall's tau : tq < tq* ■ Note that the supermodular order between U and U* 
does not necessarily yield the stochastic dominance order between C(U) and C*(U*) (i.e., 
C(U) ^ s t C*(U*) does not hold in general). For a bivariate counter-example, the interested 
reader is referred to, for instance, Caperaa et al. (1997) or Example 3.1 in Nelsen et al. (2003). 

Let us now focus on some classical families of bivariate Archimedean copulas. In Table [2j we 
obtain analytical expressions of the Kendall distribution function for Gumbel, Frank, Clayton 
and Ali-Mikhail-Haq families. 
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Copula 


9 £ 


Kendall distribution K(a, 9) 


Gumbel 
Frank 
Clayton 
Ali-Mikhail-Haq 


[l,oo) 
(-oo,oo)\{0} 
[-l,oo) \{0} 
[-1,1) 


a (l — a In a) 
«+Hl-^)m(^) 
«(l + ^(l-« 9 )) 

a -l+6M-(l-0+0a)(ln(l-0+0 a)+ln a) 


0-1 



Table 2: Kendall distribution in some bivariate Archimedean cases. 



Remark 1 Bivariate Archimedean copula can be extended to d- dimensional copulas with d > 2 
as far as the generator <j> is a d-monotone function on [0,oo) (see McNeil and Neslehovd, 2009 
for more details). For the d-dimensional Clayton copulas, the underlying dependence parameter 



must be such that 6 > 



(see Example 4-27 in Nelsen, 1999). Frank copulas can be extended 



to d-dimensional copulas for 9 > (see Example 4-24 in Nelsen, 1999). 

Note that parameter 6 governs the level of dependence amongst components of the underlying 
random vector. Indeed, it can be shown that, for all Archimedean copulas in Table [2j an increase 
of 6 yields an increase of dependence in the sense of the supermodular order, i.e., 9 < 9* =^ 
U -< S m U* (see further examples in Joe, 1997 and Wei and Hu, 2002). Then, as a consequence 
of Proposition the following comparison result holds 



< 



c(u) < sl c*(ir 



In fact, a stronger comparison result can be derived for Archimedean copulas of Table [2j as shown 
in the following remark. 

Remark 2 For copulas in Table 2, one can check that 9K Qg'^ < 0, for all a € (0,1). This 
means that, for these classical examples, the associated Kendall distributions actually increase 
with respect to the stochastic dominance order when the dependence parameter 9 increases, i.e., 

6<9* =>C(U) < st C*(U*). 

In order to illustrate this property we plot in Figure^the Kendall distribution function K(-, 9) for 
different choices of parameter 9 in the bivariate Clayton copula case and in the bivariate Gumbel 
copula case. 



2. Multivariate Value-at-Risk 

From the usual definition in the univariate setting we know that the quantile function Qx(ot) 
provides a point which accumulates a probability a to the left tail and 1 — a to the right tail. 
The univariate quantile function Qx is used in risk theory to define an univariate risk measure: 
the Value-at-Risk. This measure is defined as 

VaR a (X) = Qx(a), V a £ (0, 1). 
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Figure 1: Kendall distribution K(-,9) for different values of 9 in the Clayton copula case (left) and the 
Gumbel copula case (right). The dark full line represents the first diagonal and it corresponds to the 
comonotonic case. 

In the recent literature, an intuitive and immediate generalization of the VaR measure in the case 
of a d— dimensional loss distribution function F is represented by its a-quantile curves. More 
precisely, let F : IR^ — > [0, 1] be a partially increasing multivariate distribution function^] and 
for any a G (0,1), let L(a) = {x G : F(x) > a} be the upper a-level set of F. Tibiletti 
(1993), Embrechts and Puccetti (2006) and Nappo and Spizzichino (2009) propose to define the 
multivariate Value-at-Risk as dL(a), for a G (0, 1), where dL(a) denotes the boundary associated 
with the set L{a). 

In the following, we will consider non-negative absolutely-continuous random vectoi[^] 
X = (X\, . . . ,Xd) (with respect to Lebesgue measure A on ~R d ) with partially increasing mul- 
tivariate distribution function F and such that E(Xj) < oo, for i = 1,.. .,d. These conditions 
will be called regularity conditions. 

However extensions of our results in the case of multivariate distribution function on the entire 
space M. d or in the presence of plateau in the graph of F are possible. Starting from these 
considerations, we introduce here a multivariate generalization of the VaR measure. 

Definition 2.1 Consider a random vector X satisfying the regularity conditions. For a G (0, 1), 
we define the multidimensional Value-at-Risk at probability level a by 



4 A function F(x\, . . . , Xd) is partially increasing on \ (0, . . . , 0) if the functions of one variable g(-) = 
F(xi, . . . ,Xj-i,-,Xj+i, . . . ,Xd) are increasing. About properties of partially increasing multivariate distribution 
functions we refer the interested reader to Rossi (1973), Tibiletti (1991). 

5 We restrict ourselves to because, in our applications, components of d— dimensional vectors correspond to 
random losses and are then valued in K + . 
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Analogously 

VaR Q (X) = E[X\ F(X) = a] 
where dL(a) is the a-level set of F. 



I E[X 1 \F(X) = a] \ 
V E[X d \F(X) = a] J 



From now on, we denote by VaR^,(X), . . ., VaR^X) the components of the vector VaR a (X). 
If X is an exchangeable vector, then VaR^(X) = VaR^(X), for i,j = l,...,d. Furthermore, 



given a random variable X, since E[X \X = VaR a (X)] = Qx(a), for all a in (0, 1), Definition 2.1 
can be viewed as a natural multivariate version of the univariate one. Moreover, the proposed 
generalization of Value- at-Risk for multivariate portfolio can be seen as a more parsimonious 
and synthetic measure compared with Embrechts and Puccetti (2006) 's approach. Indeed, the 
measure proposed by Embrechts and Puccetti (2006) is represented by an infinite number of 
points (an hyperspace of dimension d — 1). This choice can be unsuitable when we face real risk 
management problems. In our proposition, instead of considering the whole hyperspace dL(a) 
corresponding to the a-level set of F, we only focus on the particular point in WL that matches 
the conditional expectation of X given that X falls in dL(a). This means that our measure is a 
real- valued vector with the same dimension as the considered portfolio of risks. The latter feature 
could be relevant on practical grounds. 



Note that, under the regularity conditions, dL(a) = {x G W± : F(~x) = a} has Lebesgue- measure 
zero in Ri (e.g., see Property 3 in Tibiletti, 1990). Then we make sense of Definition 
the limit procedure in Feller (1966), Section 3.2: 



2.1 



using 



E[Xi\F(X) = a} = lim E[X t |a<F(X) < a + h] 

h->0 



1 \ a+h . — ( 4 ) 



for i = 1, . . . , d. 

Dividing numerator and denominator in Q by h, we obtain, as h — > 

Iq x (a) x f(Xi,F(X))(x, a) dx 
E[ Xl |F(X) = a}= VaR^(X) = Qx ^ Ll± , ( 5 ) 

for i = 1, . . . , d, where K'{a) = d ^ a ^ is the Kendall distribution density function. This procedure 



gives a rigorous sense to our VaR a (X) in Definition 2.1 Remark that the existence of /(Xi,_F(x)) 
and K' in ^ is guaranteed by the regularity conditions (for further details, see Proposition 1 in 
Imlahi et al., 1999 or Proposition 4 in Chakak and Ezzerg, 2000). 

2.1. Archimedean copula case 

Surprisingly enough, the VaR introduced in Definition |2 . 1| can be computed analytically for any 
d— dimensional random vector with an Archimedean copula dependence structure. This is due to 
McNeil and Neslehova's stochastic representation of Archimedean copulas. 
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Proposition 2.1 (McNeil and Neslehovd, 2009) Let U = (U\, . . . ,Ud) be distributed according 
to a d- dimensional Archimedian copula with generator <f>, then 

(<l>(Ui),...,<KU d ))±RS, (6) 

where S = (Si, . . . , Sj) is uniformly distributed on the unit simplex |a; > | Ylt=i x k = 1 1 an d R 
is an independent non-negative scalar random variable which can be interpreted as the radial part 
of (4>(Ui), . . . , 4>(Ud)) since Ylt=i^k = 1- The random vector S follows a symmetric Dirichlet 

distribution whereas the distribution of R = Ylt=i * s directly related to the generator <f> 

through the inverse Williamson transform of (j)^ 1 . 

Recall that a (f-dimensional Archimedean copula with generator <f> is defined by C(u\, . . . ,Ud) = 
4>^ 1 (4>(u\) + • • • + 4>(ud)), for all (m, . . . ,Ud) G [0, l] d . Then, the radial part R of representation 
Q is directly related to the generator <f> and the probability integral transformation of U, that 
is, 

R = (/)(C(\J)). 

As a result, any random vector (U\, . . . , Ud) which follows an Archimedean copula with generator 
4> can be represented as a deterministic function of C(U) and an independent random vector 
S = (Si, . . . , Sd) uniformly distributed on the unit simplex, i.e., 

([/!,..., U d ) I (r 1 (Srf (C(U))) , . . . , r 1 (S d <f> (C(U)))) ; (7) 

The previous relation allows us to obtain an easily tractable expression of VaR a (X) for any 
random vector X with an Archimedean copula dependence structure. 

Corollary 2.1 Let X be a d-dimensional random vector with marginal distributions F\, . . . ,Fd- 
Assume that the dependence structure of X is given by an Archimedian copula with generator (p. 
Then, for any k = 1, . . . , d, 

VaR k a (X) = E [F- 1 (rHSkHa)))] (8) 
where Sk is a random variable with Beta(l,d— 1) distribution. 

Proof: Note that X is distributed as (F 1 _1 (C7i), . . . , F d ~ 1 (Ud)) where (U\, . . . ,Ud) follows an 
Archimedean copula with generat or <j). Then, each component k = 1, . . . , d of the multivariate risk 

can be expressed as VaR^(X) = E [^(£4) | C(U) 



measure introduced in Definition 



2.1 



a 



Moreover, from representation Q the following relation holds 

[U | C(U) = a] = (r l (5i0 (a)) , . . . , JT 1 {Sd<t> (a))) (9) 

since S and C(U) are stochastically independent. The result follows from the fact that the ran- 
dom vector (Si, . . . , Sd) follows a symmetric Dirichlet distribution. □ 

Note that, using @, the marginal distributions of U given C(U) = a can be expressed in a very 
simple way, that is, for any k = 1, . . . , d, 

F(U k < u I CCU) = a) = ( 1 - ^4^1 for < a < u < 1. (10) 

V n a )J 
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The latter relation derives from the fact that Sk, which is Beta(l, d — 1)- distributed, is such that 
Sk = 1 — V d ~ 1 where V is uniformly-distributed on (0, 1). 

The following of this section is devoted to calculate the analytical expression of VaR^ (X) in some 
Archimedean multivariate cases, using Q. 

Clayton family in dimension 2: 

As a matter of example, let us now consider the Clayton family of bivariate copulas. This family 
is interesting since it contains the counter-monotonic, the independence and the comonotonic 
copulas as particular cases. Let (X, Y) be a random vector distributed as a Clayton copula with 
parameter 9 > — 1. Then, X and Y are uniformly-distributed on (0, 1) and the joint distribution 
function Cq of (X, Y) is such that 

C e (x,y) = (m a x{x- e + y- e -l^})- 1 / , for > -1, (x, y) € [0, l] 2 . (11) 

Since X and Y are exchangeable, the two components of the multivariate VaR are identical. Using 
Corollary [2~T| Tabled gives analytical expressions for the first (equal to the second) component 



of the VaR, i.e., VaR* e (V,Y). Note that the latter can be represented as a function of the 
risk level a and the dependence parameter 9. For 9 = — 1 and 9 = oo we obtain the Frechet- 
Hoeffding lower and upper bounds: W(x, y) = max{x + y — 1,0} (counter-monotonic copula) and 
M(x, y) = mm{x, y} (comonotonic random copula) respectively. The settings 9 = and 9 = 1 
correspond to degenerate cases. For 9 = we have the independence copula H(x,y) = xy. For 



1, we obtain the copula denoted by in Nelsen (1999), where j^zjj(x,y) 



x+y—x y ' 



Copula 


9 


VaR^(X,Y) 


Clayton Cg 

Counter-monotonic W 

Independent II 
n 

s-n 

Comonotonic M 


(-l,oo) 
-1 

1 

00 


a e -a 
9-1 q»-1 

1 + Q 

2 

a-1 
In a 
a In a 
a-1 

a 



Table 3: VaR ad (X,Y) for different dependence structures. 



Interestingly, one can readily show that — q^"' 6 > and — Q g a ' e < 0, for 9 > — 1 and a G (0, 1). 
This proves that, for Clayton-distributed random couples, the components of our multivariate 
VaR are increasing functions of the risk level a and decreasing functions of the dependence para- 
meter 9. Note also that the multivariate VaR in the comonotonic case corresponds to the vector 
composed of the univariate VaR associated with each component. These properties are illustrated 
in Figure [2] where VaR^ (X, Y) is plotted as a function of the risk level a for different values of 
the parameter 8. Observe that an increase of the dependence parameter 9 tends to lower the VaR 
up to the perfect dependence case where XaR^ g (X, Y) = VaR Q (V) = a. The latter empirical 
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behaviors will be formally confirmed in next sections. 



— — 0- 1 01 9= i e=s 6 = m 

0.9 
OS 
0.7 
0.6 

0.4 
0J 
02 
0.1 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.S 0.9 

a 

Figure 2: Behavior of VaR^ e (X, Y) — VaR^ S (X, Y) with respect to risk level a for different values of 
dependence parameter 9. The random vector (X, Y) follows a Clayton copula distribution with parameter 
9. 

Ali-Mikhail-Haq in dimension 2: 



Let (X, Y) be a random vector distributed as a Ali-Mikhail-Haq copula with parameter 8 G 
[— 1, 1). In particular, the marginal distribution of X and Y are uniform. Then, the distribution 
function Cg of (X, Y) is such that 

Cg(x,y) = ^ r, for 0G [-1,1), (x, y) G [0, l] 2 . 

l-d(l-x)(\-y) 

Using Corollary [2"H Tableg gives analytical expressions for the first (equal to the second) compon- 
ent of the VaR, i.e., VaR* g(X, Y). When = we obtain the independence copula H(x, y) = xy. 



Copula 


9 


VaR^(X,Y) 


Ali-Mikhail-Haq copula Cg 
Independent n 


[-1,1) 



ln(l-0+0a)(-l+0) 


(-ln(a)+ln(l-0+0a))0 

a-1 
ln(a) 



Table 4: VaR^ g (X,Y) for a bivariate Ali-Mikhail-Haq copula. 




13 



Clayton family in dimension 3: 



We now consider a 3-dimensional vector X = (X\, X2, X%) with Clayton copula and parameter 
9 > — I (see Remark [l| and uniform marginals. In this case we give an analytical expression of 
VaR^ g (Xi, X2, X3), for % = 1,2, 3. Results are gathered in Table [HJ 



Copula 


e 


V S JV a j(X 1 ,X 2 ,X 3 ) 


Clayton Cg 
Independent II 


(-l/2,oo) 



9(-2a e-a 2e +a e +ea 2e +ea) 


(26>-l)(-l+0)(l-2a e +a 2fl ) 

1— a+\n(a) 
(Ha)f 



Table 5: VaR* g (Xi, X2, X3) for different dependence structures. 



As in the bivariate case above, one can readily show that — da a ' > and — ag"' 9 < 0, for 
9 > —1/2 and a £ (0, 1). Then, as in the bivariate case above, an increase of the dependence 
parameter 9 tends to lower the VaR up to the perfect dependence case where VaR^ d (X) = 
VaR Q (X) = a, for i = 1,2,3. These empirical behaviors will be formally confirmed in next 
sections. 



2.2. Invariance properties 

Our aim in the present section is to analyze the multivariate Value- at-Risk introduced in Defini- 
tion 2.1 in terms of classical invariance properties of risk measures (we refer the interested reader 



to Artzner et ai, 1999). The following proposition proves positive homogeneity and translation 
invariance for VaR a (X). 

Proposition 2.2 Consider a random vector X satisfying the regularity conditions. For a £ 
(0,1), X&R a (X) satisfies the following properties: 

Positive Homogeneity: V c G 



*■+> 



VaR„(cX) = cV&R n (X) 



Translation Invariance: V c E , 



VaR Q (c + X) = c + VaR Q (X) 



/ ciE[Xi|F(X) = a) 
V c d E[X d \F(X) = a 

( Cl +E[Xi\F(X) = 



a 



V c d + E[X d \F(X) = a[ 



The proof relies on the following lemma which states an invariance property for level sets. We 
refer the interested reader to Proposition 1 in Tibiletti (1993). 
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Lemma 2.1 Let the function h be such that h(x\, . . . , x d ) = (/ii(xi), . . . , h d (x d )) where hi, . . . , h d 
are non- decreasing functions, then 



( E[h 1 (X 1 )\F(X) = a[ 



VaR a (h(X)) 



V E[h d (X d )\F(X) = a] 
where F is the cumulative distribution function of X. 

Proof: Let h : — > M.+ be such that h(x\, . ■ ■ , x d ) = (hi(x{), . . . , h d (x d )) where h%, . . . , h d are 
increasing functions. Then we obtain 

9L hm {a) = {x G R d + : F h(X) (x) = a} = {x G : F (hl(Xl)> ... MXd)) (x) = q} 

= {fe(x) G R d + : F x (x) = a} = {(fci(a!i), . . . , h d (x d )) G R+ : F x (x) = a} = h{dL{a)). 

Then 

E[/ipQ) | ^(X) G dL MX) H] = E[/i(Xj) | h(X) G 5L h(x) (a)] 
= E[fc(X<)|^(X) G /i(5L x )(a)] =E[h(X i )\h~ 1 (h(X)) G 9L x (a)], 

for i = 1, . . . , d. Hence the result. □ 

2. 3. Comparison of univariate and multivariate VaR 

Note that, using a change of variable, each component of the multivariate Value-at-Risk can 
be represented as an integral transformation of the associated univariate Value-at-Risk. Let us 
denote by the marginal distribution functions of Xi for i = 1 , . . . , d and by C the copula 
associated with X. Thanks to Sklar's theorem we have F(x\, . . . , x d ) = C(Fx 1 (xi), ■ ■ ■ , Fx d (x d )) 
(see Sklar, 1959). Then the random variables Ui defined by Ui = Fx^Xi), for % = 1, . . . , d, are 
uniformly distributed and their joint distribution is equal to C. Using these notations, we get 

VaR^(X) = J^)J a VaR 7 (X 4 )/ (c/l ,C(U))(7,a)^7> ( 12 ) 

for i = l,...,d, where ffu u c(V)) 1S the density function associated with the multivariate vector 
(Ui,C(U)). The following proposition allows us to compare univariate and multivariate Value- 
at-Risk. 

Proposition 2.3 Consider a random vector X satisfying the regularity conditions. Assume that 
its multivariate distribution function F is a quasi concav^ Then, for all a G (0, 1), the following 
inequality holds 

VaRi(X) > VaRa(Xi), (13) 

fori = l,...,d. 



6 A function F is quasi concave if the upper level sets of F are convex sets. Tibiletti (1995) points out families 
of distribution functions which satisfy the property of quasi concavity. For instance all the Archimedean copulas 
are quasi concave functions (see Theorem 4.3.2 in Nelsen, 1999 for proof in dimension 2; Proposition 3 in Tibiletti, 
1995, for proof in dimension d). 
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Proof: Let a S (0,1). From the definition of the accumulated probability, it is easy to show 
dL(a) is interiorly bounded by the marginal univariate quantile functions. Moreover, recall that 
L(a) is a convex set in from the quasi concavity of F (see Section 2 in Tibiletti, 1995). Then, 
for all x = (xi, . . . , Xd) 6 dL(a), x\ > VaRo,(Xi), • • • ,Xd> V&R a (Xd) and trivially, VaR^(X) is 
greater than YaR a (Xi), for i = 1, . . . , d. Hence the result. □ 



Proposition 2.3 states that the multivariate VaR,Q,(X) is a more conservative measure than the 



vector composed with the univariate a- Value-at-Risk of marginals. Furthermore, we can prove 



that the previous lower bound in (13) is reached for comonotonic random vectors. 



Proposition 2.4 Consider a comonotonic non-negative random vector X. Then, for all a 6 
(0, 1), it holds that 

VaRi(X) = V&R a {Xi), 

fori = l,...,d. 

Proof: Let a £ (0,1)- If X = (X\, . . . , Xd) is a comonotonic non- negative random vector 
then there exist a random variable Z and d increasing functions gi, ■ ■ ■ ,gd such that X is equal 
to (gi(Z), . . . , gd{Z)) in distribution. So the set {(x±, . . . ,Xd) ■ F(x\, . . . , Xd) = a} becomes 
{(x\, . . . ,Xd) ■ m.m{g^ 1 (xi), . . . , g^ 1 (xd)} = Qz{ot)}, where Qz is the quantile function of Z. 
Then, trivially, VaR^(X) = E[ Xi \ F(X.) = a] = Qxi{o>), for i = 1, . . . , d. Hence the result. □ 



Remark 3 For bivariate independent random couple (X, Y), formula (12) becomes 



VaRi(*, y) = ^Vl / 1V0E ^ 
-ln(a)7 Q 7 

then, obviously, in this case the X-related component only depends on the marginal behavior of 
X. For further details the reader is referred to Corollary 4-3.5 in Nelsen (1999). 

2.4. Behavior of the multivariate VaR with respect to marginal distributions 

In this section we study the behavior of our risk measure with respect to a change in marginal 
distributions. Results presented below provide natural multivariate extensions of some classical 
results in the univariate setting (see, e.g., Denuit and Charpentier, 2004). 

Proposition 2.5 Let X and Y be two d-dimensional continuous random vectors satisfying the 
regularity conditions and with the same copula structure C. If Xi = Y,- t then it holds that 

VaR^(X) = VaR* Q ( Y) , for all a G (0, 1). 



The proof of the previous Proposition directly comes down from formula (12). From Proposition 



2.5 



we remark that, for a fixed copula structure C, the i-th component VaR^(X) does not depend 
on marginal distributions of the other components j with j 7^ i. 

In order to derive the next result, we use the definitions of stochastic orders presented in Section [T] 
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Proposition 2.6 Let X and Y be two d-dimensional continuous random vectors satisfying the 
regularity conditions and with the same copula structure C. If X; ^ s t Y{ then it holds that 

VaR^(X) < VaR l Q ( Y) , for all a 6 (0, 1). 



The proof comes down from formula (12) and Definition |1,1| Note that, the result in Proposition 
2.6 is consistent with the one-dimensional setting (see Section 3.3.1 in Denuit et al, 2005). Indeed, 
as in dimension one, an increase of marginals with respect to the first order stochastic dominance 
yields an increase in the corresponding components of VaR a (X). 

As a result, in an economy with several interconnected financial institutions, capital required for 
one particular institution is affected by its own marginal risk. But, for a fixed dependence struc- 
ture, the solvency capital required for this specific institution does not depend on marginal risks 
bearing by the others. Then, our multivariate VaR implies a "fair" allocation of solvency capital 
with respect to individual risk-taking behavior. In other words, individual financial institutions 
may not have to pay more for risky business activities undertook by the others. 

2. 5. Behavior of multivariate VaR with respect to the dependence structure 

In this section we study the behavior of our risk measure with respect to a variation of the 
dependence structure, with unchanged marginal distributions. 

Proposition 2.7 Let X and X* be two d-dimensional continuous random vectors satisfying the 
regularity conditions and with the same margins Fx i and Fx*, for % = 1,.. .,d, and let C (resp. 
C*) denote the copula function associated with X {resp. X*). Let Ui = Fx^Xi), U* = Fx^iX*), 
U=(U 1 ,...,Ud) and U* = (Uf, . . . ,U*). 

If [Ui\C{U) = a] < st [U*\C*(U*) = a] then VaR l Q (X) < VaR^(JT). 

Proof: Let U x = [Ui\C(U) = a] and U 2 = [U*\C*{\]*) = a]. We recall that U x < st U 2 if and 
only if E[/(t r i)] < E[/([/2)], for all non-decreasing function / such that the expectations exist 
(see Denuit et al., 2005; Proposition 3.3.14). We now choose f(u) = Qx t (u), for u G (0, 1). Then, 
we obtain 

E[QxM)\C(V) = a] < E[Q Xi ([/*)|C*(U*) = a], 
But the right-hand side of the previous inequality is equal to E[ Qx* (U*)\C*(XJ*) = a] since Xi 



and X* have the same distribution. Finally, from formula (12) we obtain VaR^(X) < VaR^(X*). 
Hence the result. □ 



We now provide an illustration of Proposition 2.7 in the case of d— dimensional Archimedean 
copulas. 

Corollary 2.2 Consider a d-dimensional random vector X, satisfying the regularity conditions, 
with marginal distributions Fx t , for i = l,...,d, and copula C. If C belongs to one of the d- 
dimensional family of Archimedean copulas introduced in Table [£| an increase of the dependence 
parameter 9 yields a decrease in each component ofV&R a (X). 
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Proof: Let Cg and Cg* be two Archimedean copulas of the same family with generator <fi@ and cj)$* 
such that 9 < 9*. Given Proposition 2.7 we have to check that the relation [U*\Cg* (TJ*) = a] < s t 



[Ui\Ce(U) = a) holds for alH = 1, . . . , d where (Ui, . . . , U d ) and (U^, . . . , U^) are distributed 
(resp.) as Cg and Cg*. However, using formula ( flO| ), we can readily prove that the previous 
relation can be restated as a decreasing condition on the ratio of generators 4>g* and cj)g, i.e., 

(j)g* 

{U*\Cg*(XJ*) = a] zi s t [Ui\Cg(U) = a] for any a G (0, 1) <^=^ — — is a decreasing function. 

(f>e 

Eventually, we have check that, for all Archimedean family introduced in Table [2l the func tion 



defined by ^ L is indeed decreasing when 9 < 9* . We immediately obtain from Proposition 2.7 
that each component of VaR a (X) is a decreasing function of 9. □ 

Then, for copulas in Table [2j the multivariate VaR is non-increasing with respect to the de- 
pendence parameter 6 (coordinate by coordinate). In particular, this means that, in the case 
of Archimedean copula, limit behaviors of dependence parameters are associated with bounds 
for our multivariate risk measure. For instance, let (X, Y) be a bivariate random vector with 
a Clayton dependence structure and fixed margins. If we denote by VaR| a ^ (X, Y) the first 
component of VaR when the dependence parameter is equal to 9, then the following comparison 
result holds for all a £ (0, 1) and all 9 G (— l,oo) 

VaRf ai00) (X, Y) < VaR^X, Y) < VaR^pT, Y). 

Note that the lower bound corresponds to comonotonic random variables, so that VaR^ a ^(X, Y) = 
VaR Q (X) = a for random variables X, Y with uniform margins (see Table [3]). The upper bound 
corresponds to counter-monotonic random variables, so that VaR^, _ 1 ^(X, Y) = for random 
variables X, Y with uniform margins, which turns to be also equal to CTE Q! (X) in that case, where 
CTE stands for the univariate Conditional- Tail- Expectation defined in Section [3] (see (|15|)). 

2.6. Behavior of multivariate VaR with respect to risk level 

In order to study the behavior of the multivariate Value-at-Risk with respect to risk level a, 
we need to introduce the positive regression dependence concept. For a bivariate random vector 
(X, Y) we mean by positive dependence that X and Y are likely to be large or to be small 
together. An excellent presentation of positive dependence concepts can be found in Chapter 2 
of the book by Joe (1997). The positive dependence concept that will be used in the sequel has 
been called positive regression dependence (PRD) by Lehmann (1966) but most of the authors 
use the term stochastically increasing (SI) (see Nelsen, 1999; Section 5.2.3). 

Definition 2.2 (Positive regression dependence) A bivariate random vector (X, Y) is said 
to admit positive regression dependence with respect to X, PRD(Y|X), if 

[Y\X = X1 ] < st [Y\X = x 2 ], Vxi<x 2 . (14) 



Clearly condition in (14) is a positive dependence notion (see Section 2.1.2 in Joe, 1997). 



From Definition 2.2, it is straightforward to derive the following result. 
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Proposition 2.8 Consider a d- dimensional random vector X, satisfying the regularity condi- 
tions, with marginal distributions Fx i; for i = l,...,d, and copula C. Let Ui = Fx i (Xi) and 
U = (Ui, . . . , Ud)- Then it holds that : 

If (Ui,C{U)) is PRD(Ui\C(U)) then VaR^(X) is a non-decreasing function of a. 



The proof of Proposition 2.8 essentially comes down from formula (12) and Definition 2.2 



Proof: We have [C/j|C(U) = a\] ^ s t [Ui\C(U) = a 2 ], Vai < a 2 - As in the proof of Proposition 

nQxM)\C(U) =ai] <E[Q Xi (Ui)\C(\J) = a 2 ). 

Then VaR^X) < VaR^ 2 (X), for any a\ < 02 which proves that VaR^,(X) is a non-decreasing 
function of a. □ 

Note that behavior of the multivariate VaR with respect to a change in the risk level does not 
depend on marginal distributions of X. 



The following result proves that assumptions of Proposition 2.8 are satisfied in the large class of 
(i-dimensional Archimedean copulas. 

Corollary 2.3 Consider a d-dimensional random vector X, satisfying the regularity conditions, 
with marginal distributions Fx t , for i = 1, . . . , d, and copula C . 

If C is a d-dimensional Archimedean copula then VaR^(X) is a non- decreasing function of a. 



Proof: From formula (10), we remark that ¥[Ui > u \ C(U) = a] is a non-decreasing function of 
a. The result then derives from Proposition 2.8 □ 



3. Multivariate Conditional- Tail- Expectation 

As well as in the univariate case, the multivariate VaR at a predetermined level a does not 
give any information about the thickness of the upper tail of the distribution function. This is 
a considerable shortcoming of the VaR measure because in practice we are not only concerned 
with the frequency of the default but also with the severity of loss in case of default. In order 
to overcome this problem, another risk measure has recently received growing attentions in the 
insurance and finance literature: the Conditional- Tail- Expectation (CTE). Following Artzner et 
al. (1999) and Dedu and Ciumara (2010), for a continuous loss distribution function F x the CTE 
at level a is defined by 

CTE a (X) =E[X\X > VaR a (X)), (15) 

where VaR a (X) is the univariate Value-at-Risk introduced above. For a comprehensive treatment 
and for references to the extensive literature on V&R a (X) and CTE Q (X) one may refer to Denuit 
et al. (2005). 

In the following, we propose a multivariate generalization of the bivariate Conditional- Tail- 
Expectation, previously introduced by Di Bernardino et al. (2011). 
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Definition 3.1 Consider a d-dimensional random vector X satisfying the regularity conditions. 
For a £ (0, 1), we define the multivariate a- Conditional- Tail- Expectation by 



CTEJX) = E\X\X£ L(a) 



Analogously, 



CTE Q (X) = E[X\F(X) > a] 



I E[Xi \Xe L{a)] 
\ E[X d \X€L(a)} 
/ E[X 1 \F(X) > a] 
\ E[X d \F(X) >a] 



From now on, we denote by CTE* (X), . . ., CTE^(X) the components of the vector CTE Q (X). 
Note that this multivariate Conditional- Tail- Expectation is a natural extension of the univariate 
one. Moreover, if X is an exchangeable vector, CTE^,(X) = CTE^(X), for i,j = 1, . . . ,d. 

From Definition 3.1 and formula ^ we write, for a G (0, 1) 



(te;,(X)- / x M^---^d) dxi ,,, dXd 

xeL(a) A(aj 



(j^x l (7) :r/ (^^(X))(^ 5 7)dx) d 7 



;«^(7)d7 



• (16) 



Formula (16) will be useful in Proposition 3.3 and Corollary 3.2 below. 



Bivariate Archimedean copula case 

We consider here a random couple (X, Y) which follows a Clayton copula distribution with 
parameter 6 > —1 as in (11). Here, X and Y are uniformly distributed. We obtain in Table [6] a 
closed- form expression for the multivariate CTE in that case. 



Copula 


9 


CTE^(X,y) 


Clayton Cq 


(-l,oo) 


1 6-l-a 2 (l+6)+2a 1+e 

2 e-i e-a{i+e)+a 1 + e 


Counter-monotonic W 


-1 


1 l-a 2 +2 In a 
4 1— o+lna 


Independent n 





1 (1-a) 2 

2 1— a+alna 


n 


1 


1 l+a 2 (2 lna^l) 


s-n 


2 (1-a) 2 


Comonotonic M 


00 


1+a 
2 



Table 6: CTE^ g(X, Y) for different copula dependence structures 

aCTEi 



Interestingly, one can readily show that 



— — > and — gf^- < 0, for 9 > -1 and a G 
(0, 1). This proves that, for Clayton-distributed random vectors, the components of our CTE are 
increasing functions of the risk level a and decreasing functions of the dependence parameter 9. 
Note also that, in the comonotonic case, our CTE corresponds to the vector composed of the 
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univariate CTE associated with each component. These properties are illustrated in Figure [3] 
where CTE* e (X,Y) is plotted clS cl function of the risk level a for different values of 9. Note 
that we obtained exactly the same feature for the multivariate VaR (see Figure [2]). The previous 
empirical behaviors will be formally investigated in next sections. 



9 = -1 6 = -0.1 6= 1 6 = S 8=oo Var_a(X) 




02 0.4 0.6 O.S 



Figure 3: Behavior of CTE* e (.X",Y) = CTE^ g (X, Y) with respect to risk level a for different values of 
dependence parameter 9. The random vector (X, Y) follows a Clayton copula distribution with parameter 
9. The horizontal line corresponds to E[X] = |. 



For analogous results in the case of Gumbel copula the interested reader is referred to Appendix. 
3.1. Invariance properties 



As in Section 2.2, we analyze here the multivariate Conditional- Tail- Expectation in terms of 
classical invariance properties of risk measures (we refer the interested reader to Artzner et al, 
1999). The following proposition proves positive homogeneity and translation invariance for 
CTE Q (X). 

Proposition 3.1 Consider a random vector X satisfying the regularity conditions. For a £ 
(0,1), CTE Q (X) satisfies the following properties: 



Positive Homogeneity: V c 6 

CTE Q (cX) = cCTE a (X) 



/ c 1 E[X 1 \F(X) > a] 
V c d E[X d \F(X) > a] 
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Translation Invariance: V c € Mi, 



CTE a (c + X) = c + CTE a {X) 



( c 1 +E[X 1 \F(X) > a) 



V c d + E[X d \F(X)>a] 



Arguments completely analogous to those in the proof of Proposition 2J2 are used to prove 
Proposition |3.1| 



Remark 4 For a = 0, using both the definition of CT~E a (X) and the definition of the a-upper 
level set L(a), we obtain 

( E[AY 



CTEq(X) 



ELX1. 



V nx d ] 

Then, as in the univariate case, the multivariate CTE is equal to the expected value of the un- 
derlying random vector for a = 0. 

3.2. Comparison with VaR and univariate CTE 

In the following, we prove a comparison result between CTE a (X) and univariate Value-at-Risk 



(see Proposition 3.2). 



Proposition 3.2 Consider a d-dimensional random vector X satisfying the regularity condi- 
tions. Assume that its multivariate distribution function F is quasi concave. Then, for any 
i = 1, . . . , d, the following inequality holds: 

CTE^(X) > V&R a (Xi), for all a G (0, 1). 



The proof is completely analogous to the proof of Proposition 2.3 



Let us now compare multivariate CTE (see Definition 3.1) with multivariate VaR (see Definition 



2.1). 



Proposition 3.3 Consider a d-dimensional random vector X satisfying the regularity condi- 
tions. If VaR^(X) is a non- decreasing function of a, then it holds that 

CTE^(X) > VaR4(X), for all a E (0, 1). 

Proof: Let us remark that, as in the univariate case, the multivariate Conditional-Tail-Expectation 
can be represented as an integral transformation of the multivariate Value-at-Risk: 



CTE J Q (X) 



1 



K(a) 



K'(7)VaR;(X)d 7 . 



(17) 



Since K' is non-negative and VaR* (X) is assumed to be non-decreasing with respect to 7, we 
obtain 
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CTE^(X) > ^ £ K'( 7 )VaRi(X)d 7 = K %« H VaR^(X) = VaR^(X), 



for a G (0,1). □ 



Analogously to Proposition 2.4 we can prove that, for comonotonic random vectors, components 



of the multivariate CTE are equal to univariate CTE of the corresponding marginals. 

Proposition 3.4 Consider a comonotonic non-negative d- -dimensional random vector X. Then, 
for all a £ (0, 1), it holds that 

CTE^(X) = CTE Q (AU 

fori = l,...,d. 

Proof: Let a G (0,1). If X = (X± , . . . , Xj) is a comonotonic non- negative random vector 
then there exist a random variable Z and d increasing functions gi, ■ ■ ■ ,gd such that X is equal 
to (gi(Z), . . . , gd(Z)) in distribution. So the set {(zi, . . . , Xd) ■ F(x\, . . . , xj) > a} becomes 
{( Xl ,...,x d ) : min^rVi),...,^ 1 ^)} > Qz(a)}. Then, trivially, CTE^(X) = E[Xi\Xi > 
Qxi(oi) ], for i = 1, . . . , d. Hence the result. □ 



Remark 5 Consider a random couple (X, Y) with standard uniform marginals and independent 
copula. For all a G (0, 1), it holds that 

CTEl(X,Y) > CTE Q (X), 

CTE 2 a {X,Y) > CTE Q (y). 



Indeed CTE^(X, Y) - CTE Q (X) = -\ a ( ~ 2 _"^"J""^ lna) > 0, for all a £ (0, 1). Then in this 
particular case, coordinate by coordinate, CTE a (X, Y) is a more conservative measure than the 
usual univariate CTE of marginals. 



3.3. Behavior of multivariate CTE with respect to marginal distributions 

In this section we study the behavior of the multivariate Conditional- Tail- Expectation in Defi- 



nition 3.1 with respect to a variation of marginals. Results presented below provide a natural 
multivariate extension of results in the univariate setting (see, e.g., Denuit and Charpentier, 
2004). 



Analogously to Proposition 2.5, we can state the following result. 



Proposition 3.5 Let X and Y be two d-dimensional random vector satisfying the regularity 
conditions and with the same copula structure C. If Xi = Yi, then it holds that 

CTE^(X) = CTE* Q (Y), for all a £ (0, 1). 
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Proof: Let Fi(x) = Fi(xi, . . . ,x d ) = C(F Xl (xi),..., F Xd (x d )). We recall that 



CTE^(X) = / Xl /x(^ ---^d) ^ . . . ^ 
i xei(c) A(aJ 

We remark that K(a) depends solely on the copula structure (for further details see Theorem 9 
in Belzunce et al, 2007). Then, for a fixed copula C, K{a) is invariant to a change of marginals 
(see Section[T]). Let u := (u\, . . . ,u d ). Using a change of variables, we obtain: 

CTEi(X) = f f Q Xt{Ul )^^C(u)du ) 

\Ju:C(u)>a Of-d d J K(a) 

for a G (0, 1). Hence the result. □ 

In the following we analyze how the CTE a (X) behaves when the marginal random variables 
increase with respect to some particular stochastic orders (see, e.g., Section 3.3 in Denuit et al., 
2005). Using Definitions 1.1 1.2 and Lemma |l.l| in Section[T]we can state the following result. 



Proposition 3.6 Let X and Y be two d- dimensional random vector satisfying the regularity 
conditions and with the same copula structure C. If Xi <£, Y{ then it holds that 

CTE^(X) < CTEi( Y), for all a G (0, 1). 
Proof: Let us consider the function A defined by 
a i — y A(a) = K{a) (CTE l a (Y) - CTE^(X)) 

= ( / (Qy ! K)-^W) 7r ^ r C(u)du ] =L-, 

\Ju:C{u)>a Ol...dd J K[a) 

for a G (0, 1). Since Xi < D Fj yields ELY;] < E[Yi] then A(0) = E[3^] - E[Xi] > and obviously 
lim^i A(i) = 0. Furthermore, from the definition of the dangerousness order relation, we observe 
that there exists some real number c such that Qxi(t) > QYi(t)i for all t £ (0,Fxi(c)) and 
Qxi(t) < QYi(t), for all t G (Fxtic),!). So A first increases on (OjFx^c)) and then decreases 
on (Fx i (c),l). It thus remains non-negative and ensures that CTE^,(X) < CTE^(Y), for all 
a G (0, 1). Hence the result. □ 



From Proposition 3.6 if we analyze the marginal behavior without varying the dependence struc- 
ture we have that CTE a (X) associated with a more dangerous marginal distribution will be 
greater. This means that this measure induces a greater coverage in the case of more dangerous 
risk (in the sense of the dangerousness order). Note that if Xi Yi for all i = l,...,d, we 
obtain CTE a (X) < CTE a (Y), for all a G (0, 1). 

In the following example we provide an illustration of Propositions 3.5 and 3.6 

Example 1 We consider a bivariate Clayton copula with parameter 1 and five different bivariate 
random vectors (X,Yi), for i = 1, . . . , 5. Let X ~ Exp(l), Y\ ~ Exp(2), Y 2 ~ Burr(2,l), 
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Y% ~ Exp(l), I4 ~ Frechet(4) and I5 ~ Burr(4, 1). We assume the same copula structure for 
all vectors. We calculate CTE a (X, Yj), for i = 1, . . . , 5. The results are gathered in Table 
[7j As proved in Proposition |3.5| we obtain an invariant property on the first coordinate of all 
CTE a {X,Yi), for i = 1,...,5. 

Remark that CTE a (X, I3) is a case of an exchangeable continuous random vector then we 
have CTE^(X,y 3 ) = CTE 2 a (X,Y 3 ). Furthermore, as Q X (a) = 2Q n (o), then CTE^X^) = 
2CTE 2 a (X,Y 1 ), for a e (0,1) (see Table Q. 

Since Y% ■<£> I5 ■<£> I4 <d Y2, then, for any level a G (0, 1), 

CTE^X,^) < CTE 2 a (X,Y 5 ) < CTE^X,^) < CTE 2 a (X,Y 2 ). 

Notice that the relation <£, is not transitive, since if F and G cross once, and G and H cross 
once, then F and H may cross twice. However, in the simple cases considered in Example [TJ 
this situation does not happen and therefore the transitive relation is valid. Conversely Y^,Y^ 
and 13,15 are not ordered in dangerousness sense, and also (CTE^(X, I3), CTE^(X, I4)) and 
(CTE^(X, I3), CTE^(X, I4)) are not ordered for any level a 6 (0,1). As before we can also 
verify that CTE Q (X,Yi) < CTE Q (X,y 3 ) < CTE Q (X,F 2 ), for any level a G (0,1). 



a 


CTE a (X,li) 


CTE Q (X,y 2 ) 


CTE Q (X,y 3 ) 


CTE a (X,y 4 ) 


CTE a (X,Y 5 ) 


0.10 


(1.188, 0.594) 


(1.188, 1.838) 


(1.188, 1.188) 


(1.188, 1.315) 


(1.188, 1.229) 


0.24 


(1.449, 0.724) 


(1.449, 2.218) 


(1.449, 1.449) 


(1.449, 1.431) 


(1.449, 1.366) 


0.38 


(1.727, 0.864) 


(1.727, 2.661) 


(1.727, 1.727) 


(1.727, 1.555) 


(1.727, 1.506) 


0.52 


(2.049, 1.025) 


(2.049, 3.235) 


(2.049, 2.049) 


(2.049, 1.704) 


(2.049, 1.667) 


0.66 


(2.454, 1.227) 


(2.454, 4.074) 


(2.454, 2.454) 


(2.454, 1.902) 


(2.454, 1.876) 


0.80 


(3.039, 1.519) 


(3.039, 5.591) 


(3.039, 3.039) 


(3.039, 2.219) 


(3.039, 2.202) 


0.90 


(3.768, 1.884) 


(3.768, 8.175) 


(3.768, 3.768) 


(3.768, 2.675) 


(3.768, 2.665) 


0.99 


(6.102, 3.059) 


(6.102, 26.59) 


(6.102, 6.102) 


(6.102, 4.813) 


(6.102, 4.811) 



Table 7: GTE a (X, Yi), for i = 1,...,5, with the same copula Clayton copula with parameter 1, 
X - Exp(l) and Y x - Exp(2); Y 2 ~ Burr(2, 1); Y 3 ~ Exp(l); Y 4 - Frechet(4); Y 5 - Burr(4, 1). 



3.4. Behavior of multivariate CTE with respect to the dependence structure 



In analogy with Section 2.5, in this section we study the behavior of our CTE with respect to a 



variation of the dependence structure, with unchanged marginal distributions. 

Proposition 3.7 Let X and X* be two d-dimensional continuous random vectors satisfying the 
regularity conditions and with the same margins F Xi and Fx*, for i = 1, . . . ,d, and let C (resp. 
C*) denote the copula function associated with X (resp. X*). Let Ui = Fx t (Xi), U* = Fx *(X*), 
U=(U 1 ,...,U d ) and U* = (&,... ,U*). 

If [Ui\C{U) > a] ± st [U*\C*(IT) > a] then CTE^X) < CTE^(X*). 
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Proof: We recall that U\ U% if and only if E[/(J7i)] < E[/(C/2)], for all non-decreasing 
function /, such that the expectations exist (see Denuit et at, 2005; Proposition 3.3.14). We now 
choose f(u) = Qxi(u), for u £ (0, 1). Then we obtain 

E[Q X M)\C(V) >a}< E[Q Xi (U:)\C*(V*) > a], 
But the right-hand side of the previous inequality is equal to E[ Qx* (c7*)|C*(U*) > a] since Xj 



and X* have the same distribution. Finally, from formula ([16]) we obtain CTE^(X) < CTE* a (X*). 
Hence the result. □ 



As in Corollary 2.2 we now provide an illustration of Proposition 3.7 in the case of bivariate 



Archimedean copulas. 

Corollary 3.1 Consider a bivariate random vector X, satisfying the regularity conditions, with 
marginal distributions Fx i; for i = 1,2, and copula C. If C belongs to one of the bivariate family 
of Archimedean copulas introduced in Table [1] an increase of the dependence parameter 9 yields 
a decrease in each component of CT~E a (X). 

Proof: In the bivariate Archimedean case, the joint distribution of (U, C(U, V)) can be obtained 
analytically by using a change of variable transformation^] from (U, V) to (U, C{U, V)): 

Ma) Mu) 

F {um , V )){u,a) = a- ^ + < a < u < 1. (18) 



Then, thanks to formula (18), we can obtain the joint survival probability 



h(u, a) := P[U > u , C(U, V) > a] = 1 - u + for < a < u < 1. (19) 

(f)'(a) 

Note that K(a) = ¥[C(U, V) > a] = h(a,a). Let Cq and Cq* be two bivariate Archimedean 
copulas of the same family with generator (ftg and (/)$* such that 9 < 9* . Given Proposition 
3.7 and by exchangeability, we only have to check that the relation [U*\Cg*(U* , V*) > a] ^ s t 



[U\Ce(U, V) > a] hold where (U, V) and (U*,V*) are distributed according to (resp.) Cq and 



Cq*. Given formula (19), the previous relation can be restated as 



h*(u,a) h(u,a) 

.J < -rf^4, for < a < u < 1, 20 
n*[a,a) n(a,a) 



where, from (19), h(u,a) = l — u + 4>Q(u)/4>' e (a) and h*(u,a) = l — u + 4>Q*(u)/(j)' et (a). Eventually, 



we have checked that, for all Archimedean family introduced in Table p] relation (20) is indeed 
satisfied when 9 < 9*. We then immediately obtain from Proposition |3 . 7| that each component of 
CTE Q (X) is a decreasing function of 

Then, for copulas in Table [2] the multivariate CTE is non-increasing with respect to the depend- 
ence parameter 9 (coordinate by coordinate). In particular, this means that limit behaviors of 



7 In the book by Nelsen (1999) (Corollary 4.3.5), a geometrical argument is used instead to obtain the distribution 
function of (U, C(U, V)). 
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dependence parameters are associated with bounds for our multivariate risk measure in the case 
of Archimedean copula. For instance, if we denote by CTE^ a ^ (X, Y) the first component of the 
bivariate CTE for a vector (X, Y) with a Gumbel dependence structure with parameter 8, we get 
the following comparison result: 

CTE; Qi+oo) (x,y) < CTEf Qie) (x,y) < CTE[ Qil) (x,n 

for all a G (0,1) and all 6 6 (1, oo) (see Figure [6] in Appendix). Note that the lower bound 
corresponds to comonotonic random variables, so that CTE* Q +OQ ^(X, Y) = for random 
variables X, Y with uniform margins (see Table [6]) . The upper bound corresponds to counter- 
monotonic random variables, so that CTEf i\(X, Y) = \ 1 7 ct + ?\ lna for random variables X, Y 
with uniform margins (see Table [6] and Figure [6] in Appendix) . 

3.5. Behavior of multivariate CTE with respect to risk level 



As in Section 2.6 we study the behavior of the multivariate Value- at-Risk with respect to risk 
level a. As for the VaR, the relationship between CTE a (X) and the level a will be connected to 
some positive dependence concepts. 

Corollary 3.2 Consider a d- dimensional random vector X satisfying assumptions of Proposition 



3.3, then CTE^,(X) is a non- decreasing function of risk level a. 



Proof: Let us consider the i-th coordinate CTE^(X). From (17) we have 



^CTE^(X) = [CTE^(X) - VaP^(X)] . 

da K(a) 



Using Proposition 3.3 the latter expression is non- negative for any level a£ (0,1). □ 



As a result, if VaR^(X) is a non-decreasing function of a, then CTE^,(X) is in turn a non- 



decreasing function of a. Then, thanks to Corollary 2.3, coordinates of the proposed multivariate 
Conditional- Tail- Expectation are increasing with respect to risk level for any Archimedean de- 
pendence structure. 

Corollary 3.3 Consider a d-dimensional random vector X, satisfying the regularity conditions, 
with marginal distributions Fx il for i = 1, . . . , d, and copula C . 

If C is a d-dimensional Archimedean copula then CTE^(X) is a non- decreasing function of a. 



Proofs of Proposition 3.3 and Corollary 3.2 use direct calculations. We get below an alternative 



approach based to the positive dependence concepts (see Joe, 1997). 

Definition 3.2 (Right-tail increasing) A bivariate random vector (X, Y) is right-tail increas- 
ing with respect to X, RTI(y|A), if 

[Y\X > xi] < st [Y\X > x 2 ], Vxi < x 2 . (21) 
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Remark that in (21), Y is more likely to take large values if X increases. 



Definition 3.3 (Right conditional increasing) A bivariate random vector (X, Y) is right 
conditional increasing with respect to X, RCI{Y\X), if 

[Y\X = x] ± st [Y\X > x], Vx. 

For more details on these positive dependence concepts, the reader is referred to Belzunce et al. 
(2007). In particular in the case of absolutely continuous random vector (X,Y) it holds that 

PRB(Y\X) KTL(Y\X) <=^ RCI(Y|X), 



(see Lemma 15 in Belzunce et al., 2007). Then, using Definitions 3.2 3.3 we can alternatively 



prove Proposition 3.3 and Corollary |3. 2 



Remark 6 We note that, if (X, Y) has one of the positive dependence properties mentioned 
above, then all random vectors with the same copula than (X, Y) have also the same property (see 
Nelsen, 1999; Corollary 5.2.11). For instance, ifPKD(Y\X) then also PKD(a 2 (Y)\a 1 (X)), for 
any increasing functions a\ and a 2 - 

In the univariate setting the Conditional- Tail- Expectation contains a safety loading i.e., 
CTE Q (X) > E[X], Va G (0,1) (see Section 2.4.3.3 in Denuit et al, 2005). The safety loading 
should cover the fluctuations of loss experience. Corollary |3.4| below provides a similar property 
also for our multivariate CTE. 



Corollary 3.4 Under assumptions of Corollary 3.2, for all a G (0, 1), it holds that 



CTE^(X) >E[Xi], 



Remark 7 To summarize, in the class of d- dimensional Archimedean copulas it holds that 

CTEi(X) > XaK(X) > V&R a (Xi), 

for a G (0,1). Moreover, CTE^(X, Y) and VaR„(X, Y) are non- decreasing functions of a. Let 
us stress that this non- decreasing property only depends on the dependence structure and not on 
the marginal distribution functions. 

In Figure [4] and Figure [5j we provide illustrations of Remark [7] for two particular bivariate 
Archimedean families: Frank and Ali-Mikhail-Haq copulas. 
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Figure 5: Ali-Mikhail-Haq copula with standard uniform marginals, parameter 9 = —0.7 (left), parameter 
9 = 0.99 (right). 



Conclusion and perspectives 

In this paper, we proposed an extension of the classical Value-at-Risk and Conditional- Tail- 
Expectation risk measures for continuous random vectors. As in the Embrechts and Puccetti 
(2006) 's approach, the introduced risk measures are based on multivariate generalization of 
quantiles but they are able to quantify risks in a much more parsimonious and synthetic way: the 
risk of a d-dimensional portfolio is evaluated by a point in The proposed multivariate risk 
measures may be useful for some applications where risks are heterogeneous in nature or because 
they cannot be diversify away by an aggregation procedure. 

We analyzed our multivariate risk measures in several directions. Interestingly, we showed that 
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many properties satisfied by the univariate VaR and CTE can be translated to our proposed 
multivariate versions under some conditions. In particular, the proposed VaR and CTE both 
satisfy the positive homogeneity and the translation invariance property which are parts of the 
classical axiomatic properties of Artzner et al. (1999). Thanks to the theory of stochastic orders, 
we also analyzed the effect of some risk perturbations on these measures. In the same vein as for 
the univariate VaR and CTE, we proved that an increase of marginal risks yield an increase of our 
multivariate VaR and CTE. We also gave the condition under which an increase of the risk level 
tends to increase components of the multivariate extensions and we show that these conditions are 
satisfied for (i-dimensional Archimedean copulas. We also study the effect of dependence between 
risks on individual contribution to the multivariate risk measure and we prove that for different 
families of Archimedean copulas, an increase of the dependence parameter tends to lower the 
components of the VaR and of the CTE. At the extreme case where risks are perfectly dependent 
or comonotonic, our multivariate risk measures is equal to the vector composed of univariate risk 
measures associated with each component. This feature is in line with the observation made by 
Zhou (2010): "When regulating a system consisting of similar institutions, or in other words, the 
system is highly interconnected, considering a micro-prudential regulation can be sufficient for 
reducing the overall systemic risk." (Zhou, 2010). 

Due to the fact that the Kendall distribution is not known analytically for elliptical random 
vectors, it is still an open question whether components of our proposed measures are increasing 
with respect to the risk level for such dependence structures. However, numerical experiments 
in the case of Gaussian copulas support this hypothesis. More generally, the extension of the 



McNeil and Neslehova's representation (see Proposition 2.1 ) for a generic copula C and the study 
of the behavior of distribution [C/|C(U) = a], with respect to a, are potential improvements to 
this paper that will be investigated in future work. 

In a future perspective, it should also be interesting to discuss the extensions of our risk measures 
to the case of discrete distribution functions, using "discrete level sets" as multivariate definitions 
of quantiles. For further details the reader is referred, for instance, to Laurent (2003). Another 
subject of future research should be to compare our multivariate Conditional-Tail-Expectation and 
Value-at-Risk with existing multivariate generalizations of these measures, both theoretically and 
experimentally. 
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Appendix: Gumbel copula case 

We consider Gumbel bivariate family of copulas: 

C e (x,y) = e-((- l ^ B+ (- l ^ 6 )\ 

for 9 G [l,oo), (x,y) £ [0, l] 2 (e.g. see Section 3.3.1 in Nelsen, 1999) and X, Y standard uniform 
marginals. For 9 = 1 we have the independent copula C%(x,y) = H(x, y) = xy; for 9 = 00 the 
Frechet bound M(x,y) = min{x, y} (comonotonic random variables). In this case we get 



CTE a (X,Y) 



t(a,9) 
t(a,9) 



with 



t(a, 



1 \na{9+e 2 -a 2 e -a 2 8 2 )+2a? (- In a)~ 2 + 2 WM( § (9, | 6+ 1 , - In a)+2 a 2 (In a) 2 +2 a 2 (In a) 2 



2 (l+6»)lna(0-a0+a In a) ' 

where WM is the Whittaker function WM(/i, z) = exp~5 z 2:5 + "Hpy(| + v - fi, 1 + 2v; z) and Hpy is 
the Kummer's confluent hypergeometric function. 

Remark that X and Y are exchangeable, then obviously CTE* e (X, Y) = CTE 2 e (X, Y). Furthermore 

— Q a a ' B > and — Q e a ' B < 0, for 9 > 1 and a G (0, 1). In analogy with Figure 3 in Figure 6 we 
propose a graphical illustration of CTE^ e (X, Y), in Gumbel copula case, for different values of 
the dependence parameter 9 and with respect to risk level a. 



■6=1 



15 B = QD 



■ Var_«(X) 




Figure 6: Behavior of CTE* F) = CTE^ e (X,y) with respect to risk level a for different values of 
dependence parameter 9. The random vector {X, Y) follows a Gumbel copula distribution with parameter 
9. The horizontal line corresponds to E[X] = ~. 
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